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Abstract. The changes of the spin depolarization length in zinc-blende semicon- 
ductors when an external component of correlated noise is added to a static driving 
electric field are analyzed for different values of field strength, noise amplitude and 
correlation time. Electron dynamics is simulated by a Monte Carlo procedure which 
keeps into account all the possible scattering phenomena of the hot electrons in the 
medium and includes the evolution of spin polarization. Spin depolarization is studied 
by examinating the decay of the initial spin polarization of the conduction electrons 
through the D'yakonov-Perel process, the only relevant relaxation mechanism in III-V 
crystals. Our results show that, for electric field amplitude lower than the Gunn field, 
the dephasing length shortens with the increasing of the noise intensity. Moreover, a 
nonmonotonic behavior of spin depolarization length with the noise correlation time 
is found, characterized by a maximum variation for values of noise correlation time 
comparable with the dephasing time. Instead, in high field conditions, we find that, 
critically depending on the noise correlation time, external fluctuations can positively 
affect the relaxation length. The influence of the inclusion of the electron-electron 
scattering mechanism is also shown and discussed. 
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Possible utilization of the electron spin as an information carrier in electronic devices 
is an engaging challenge for future spin-based electronics (spintronics). In particular, 
in semiconductor spin-based devices, information is encoded in the electron spin state, 
transferred by the conduction electrons, and finally detected [Il[2l|3]. However, a disad- 
vantage of the use of spin degree of freedom is the fact that spin polarization is a quantity 
that is not preserved, as each initial non-equilibrium magnetization decays both in time 
and distance during the transport. Therefore, the spin depolarization length could be 
not long enough to permit a reliable manipulation, control and detection of information. 

In spintronic devices, the information stored in a system of polarized electron 
spins, can be transferred as attached to mobile carriers by applying an external electric 
field P, [21 El m El |6] . Because of increasing miniaturization, to avoid too much intense 
electric fields, which could lead the system to exhibit a strongly nonlinear physical be- 
havior, applied voltages are very low. Since, low voltages are more subjected to the 
background noise, it becomes essential to understand the influence of fluctuations of the 
electric field on the spin depolarization process. 

The presence of noise is generally considered a disturbance, especially in determin- 
ing the efficiency of a device, since strong fluctuations can affect its performance. The 
existence of fluctuations, for example, can limit the lifetime of the information stored 
in a memory cell, bother the opening (or closure) of random logic gates and cause the 
enlargement of the distribution of arrival times of signals in transmission lines. In quan- 
tum computation, a fundamental problem is the destruction of entangled states of qubits 
by interaction with the environment. This event is characterized by loss of coherence, 
which is not suitable for the design of quantum computers [7]. 

In the last decade, however, an increasing interest has been directed towards pos- 
sible constructive aspects of noise in the dynamical response of non-linear systems 
[H m [10]. The effect of an external source of noise on the electron transport in 
GaAs crystals in the presence of static and/or periodic electric fields has been studied 
[TT| [T2| [T3l [H]. Furthermore, theoretical works which discuss the way to improve the 
ultra-fast magnetization dynamics of magnetic spin systems by including random fields 
have been recently published |T5l |16l [T7] . Nevertheless, to the best of our knowledge, 
the investigation of the role of noise on the electron spin dynamics in semiconductors is 
still missing. 

In this paper, we investigate the spin relaxation process in low-doped n-type GaAs 
crystals driven by a randomly fiuctuating electric field. The electron transport is simu- 
lated by a Monte Carlo procedure which takes into account all the possible scattering 
phenomena of hot electrons in the medium, and includes the evolution of the spin po- 
larization vector. The effects caused by the addition of an external source of correlated 
noise are investigated by analyzing the change of the spin depolarization length with 
respect to the value obtained in absence of noise. Our findings show that the presence 
of a random contribution can affect the value of the decoherence length, and that the 
variation is maximum for values of the noise correlation time comparable with the char- 
acteristic time of the spin relaxation process. Morever, noise induced effects are slightly 
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reduced by the inclusion of the electron-electron Coulomb interaction. 

The spin-orbit interaction couples the spin of conduction electrons to the electron 
momentum, which is randomized by scattering with impurities and phonons. The spin- 
orbit coupling gives rise to a spin precession, while momentum scattering makes this 
precession randomly fluctuating, both in magnitude and orientation [18] . 

For electrons delocalized and under nondegenerate regime, the D'yakonov-Perel 
(DP) mechanism [191 EQ] is the only relevant relaxation process in n-type III-V 
semiconductors [3l [211 122] . In a semiclassical formalism, the term of the single electron 
Hamiltonian, which accounts for the spin-orbit interaction, can be written as 

Hso = \^-^- (1) 

It represents the energy of electron spins precessing around an effective magnetic field 
(B = MI/hbq) with angular frequency fi, which depends on the orientation of the 
electron momentum vector with respect to the crystal axes. Near the bottom of each 
valley, the precession vector can be written as 



Or = ^[Ukl - kl)± + ky{},l - fc^y + - A:J)z] (2) 
in the F- valley 



Ol = ^ [{ky - kf)y. + {k^ - kx)y + {k^ - ky)i] (3) 

in the L- valleys located along the [111] direction in the crystallographic axes, and 

nx = Px[-kyy + k,i] (4) 

in the X- valleys located along the [100] direction [21]. In equations ([2D-(|1|), ki {i = x, y, z) 
are the components of the electron wave vector. /3r, /3l and Px are the spin-orbit 
coupling coefficients. Here, we assume /32,=0.26 eV/A-2/h and /3x=0.059 eV/A-2/h, as 
recently theoretically estimated in Ref. [25], while /3r is calculated as in Ref. [2^. 

Since the quantum-mechanical description of the electron spin evolution is 
equivalent to that of a classical momentum S experiencing the magnetic field B, we 
describe the spin dynamics by the classical equation of precession motion 

^ = O X S. (5) 

dt ^ ' 

The DP mechanism acts between two scattering events and reorients the direction of 
the precession axis and of the effective magnetic field B in random and trajectory- 
dependent way. This effect leads the spin precession frequencies fi and their directions 
to vary from place to place within the electron spin ensemble. This spatial variation is 
called inhomogeneous broadening 



When a semiconductor is in contact with a spin polarization source at x = 0, in 
the absence of any magnetic field (B = 0) and with a static electric field directed along 
the x-direction (F = -Fx), the drift-diffusion model predicts an exponential decay of 
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the spin polarization S{x) = 5*0 exp(— x/L), in which L is the electric-field-dependent 
spin depolarization length [30]. The electron spin dephasing length can be seen as the 
distance that the electrons move with average drift velocity Vd within the spin lifetime 
r i30j, 

L = VdT = fiFr = ^, (6) 

m 

where n = qTp/m* is the electron mobility, q the elementary charge, Tp the momentum 
relaxation time and m* the electron effective mass. 

The inhomogeneous broadening, quantified by the average squared precession 
frequency <| f2(k) p>, together with the correlation time of the random angular 
diffusion of spin vector Tc are the relevant variables in the D'yakonov-Perel's formula 

T-'=<\n{k) |2>r,. (7) 

Since, in our case, we include the electron-electron interaction mechanism, one needs to 
distinguish between the momentum relaxation time Tp and the momentum redistribution 
time Tp, which is practically equal to Tc. This distinction is necessary because, although 
electron-electron scattering contributes to momentum redistribution, it does not directly 
lead to momentum relaxation [311 132] . 

By using equations ([6]) and ([7]), the spin depolarization length can be expressed in 
the form, 

m* <\ n{k) |2>r;' ^ ^ 

In our simulations the semiconductor bulk is driven by a fluctuating electric fleld 

F{t) = Fo + v{t) (9) 

where Fo is the amplitude of the deterministic part and ri{t) is a random term, modeled 
by a stochastic process. Here, ri{t) is modeled as an Ornstein-Uhlenbeck (OU) process, 
which obeys the stochastic differential equation [33] : 



"W + ./H^^W (10) 



dt td V td 

where td and D are, respectively, the correlation time and the intensity of the noise 
described by the OU process with autocorrelation function 

{'rj{t)r]{t')) = i^exp(— |t — t'l/r^i). ^{t) is a Gaussian white noise with zero mean 
< ^{t) >= 0, and autocorrelation function {C(t)^(t')) = S{t — t'). Within the framework 
of the Ito's calculus, the general solution of the equation ( ITOj) leads to the following 
expression for the stochastic evolution of the amplitude of the electric field 



I2D /■* t-t' 

F{t) = Fo + T]{0)e-'/^^ + J— / e-^dW{t'), (11) 



where the initial condition is ri{0) = 0, and W{t) is the Wiener process [33] . 

In a practical system, r](t) could be generated by a RC circuit driven by a source 
of Gaussian white noise, with correlation time td = {RC)~^ (see equation (JTOl) ). The 
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Gaussian white noise can be generated either by the Zener breakdown phenomenon in a 
diode or in a inversely polarized base-collector junction of a BJT, either by amplifying 
the thermal noise in a resistor |3l]. The correlation time td is tunable by using a 
diode (varicap) with a volt age- dependent variable capacitance; the noise intensity D 
can be chosen, for example, by suitably amplifying the noise produced through the 
Zener stochastic process. 
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Figure 1. Spin polarization {Sx) as a function of the distance from the injection 
point (.T = 0) obtained by applying a fluctuating field characterized by a deterministic 
component Fq and a random component with standard deviation D^l"^ ^ for different 
values of noise correlation time tu'- 10~^to, 10~^to, IO^tq, and in absence of noise, (a) 
Fq = 1 kV/cm, = 0.6 kV/cm and (b) Fq = 6 kV/cm, D^/^ = 3.6 kV/cm 



The Monte Carlo code used here follows the procedure described in Ref. [28] • The spin 
polarization vector has been incorporated into the algorithm as an additional quantity 
and calculated for each free carrier. Furthermore, the electron-electron scattering is 
included, as extensively described in Ref. [29]. All simulations are performed in a GaAs 
crystal with a free electrons concentration equal to 10^'^ cm~^ and lattice temperature 
equal to 77 K. The temporal step is 10 fs and an ensemble of 5 ■ 10^ electrons is used 
to collect spin statistics. Moreover, we assume that all donors are ionized and that the 
free electron concentration is equal to the doping concentration. All physical quantities 
of interest are calculated after a transient time of typically 10^ time steps, long enough 
to achieve the steady-state transport regime. The spin relaxation simulation starts with 
all electrons of the ensemble initially polarized ((S) = 1) along x-axis of the crystal, at 
the injection plane [xq = 0) [5l|3]. The spin lifetime r and the spin depolarization length 
L are calculated by extracting, respectively, the time and the distance from the injection 
plane of the center of mass of the electron ensemble, corresponding to a reduction of 
the initial spin polarization by a factor 1/e. 

In Figure [H we show the electron spin average polarization {^S^ as a function of the 
distance traveled by the center of mass of the electron cloud from the injection point, 
in the presence of a fluctuating field characterized by a deterministic component with 
amplitude Fq and a random component with standard deviation D^/^, for three different 
values of noise correlation time td'- 10~^ro, 10~^ro, IO^Tq and in absence of noise; tq is 
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Figure 2. (a) Ratio between tlie spin depolarization length L in the presence of noise 
and Lq, obtained in absence of noise, as a function of the ratio between the noise 
correlation time tjj and the spin relaxation time in absence of noise tq, at different 
values of noise intensity D. (b) Electric field amplitude F{t) and (c) squared precession 
frequency | fl(k) p as a function of time t. Fq = 1 kV/cm. 



the spin relaxation time obtained in absence of noise at the same value of Fq. In panel 
(a): Fo = 1 kV/cm and D^/^ ^ q q kV/cm; in panel (b): Fq = 6 kV/cm and D^/^ ^ 3 5 
kV/cm. 

When td « tq, the spin dephasing process is not affected by the fluctuations of 
the electric field, which have a negligible memory (td) with respect to the characteris- 
tic time To of the system. In particular, under low electric field condition [panel (a)], 
the spin relaxation process is significantly influenced by the fluctuating field only for 
values of noise correlation time comparable with tq, while the process becomes quasi- 
deterministic when Td » tq. Differently, under high electric field condition [panel (b)] 
when Td ~ tq and td » tq, a slow down of the spin relaxation process is observed. 

With the aim of investigating the effects of the correlated noise source on the spin 
depolarization process, we performed 100 different realizations and evaluated both av- 
erage values and error bars of the extracted spin depolarization lengths. 

In panel (a) of Figure O we show the ratio between the spin depolarization length 
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L in the presence of noise and Lq, obtained in absence of noise, as a function of the ratio 
between the noise correlation time r/) and tq, at different values of noise intensity D and 
with Fq = 1 kV/cm. For these values of parameters, Lq = 32.6 /im and tq = 0.16 ns. 
The addition of a source of correlated fluctuations, characterized by IO'^tq < td < tq, 
reduces the values of the spin depolarization length L up to 25%. In particular, L/Lq is 
a nonmonotonic function of Td/tq which exhibits a minimum for Td/tq ^ 0.1. For both 
Tf) « Tq and Tf) » Tq, the values of L coincide with those of Lq. The presence of 
the minimum, which becomes deeper with the increasing of the noise amplitude, can be 
explained by analyzing the temporal evolution of the quantities related to the electron 
transport and to the spin relaxation process. In panels (b) and (c) of Figure Ej we report 
the electric field amplitude F{t) and the squared precession frequency | 0(k) p, respec- 
tively, as a function of time t. For very low values of to/tq, \ n(k) ^ symmetrically 
fluctuates around its average value, corresponding to that obtained in absence of noise. 
By increasing the value of td, the effective electric field felt by electrons, within a time 
window comparable with the spin relaxation time, becomes very different from the value 
Fq [panel (b)]. As a consequence, the temporal evolution of | r2(k) p shows an evident 
asymmetry in the same temporal window [panel (c)]. Because of the proportionality 
between the electron momentum and the electric field F{t), the equation (|2]) leads 
to a quadratic relation between F{t) and | 0(k) p on the k'^{ky — klY term and at 
fourth power on the other two terms. Hence, the values of F greater than Fq give rise 
to values of | f2(k) p much greater than those obtained for F < Fq. So, in accordance 
with equation (^^, the asymmetry of | r2(k) p is responsible for the observed reduction 
of spin lifetime. By further increasing the value of T£,, the random fluctuating term r]{t) 
of the electric field tends to its initial value //(O) = [see equation (ITOll]. and F{t) — )■ Fq. 
Therefore, the behavior of system becomes quasi-deterministic and the spin dephasing 
length L approaches its deterministic value Lq. 
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Figure 3. (a) Ratio between tlie spin depolarization length L in the presence of 
noise and Lq, obtained in absence of noise, as a function of the ratio between the noise 
correlation time td and the spin relaxation time in absence of noise tq and (b) electron 
occupation percentage in L-valleys 77^ as a function of td/tq, at different values of noise 
intensity D. Fq = Q kV/cm. 
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In Figured (a), we show the ratio L/ as a function of T^jr^ for Fq = 6 kV/cm, 
at different values of noise intensity D. For these values of parameters, Lq = 279 nm 
and To = 1.13 ps. In the presence of a driving electric field greater than the necessary 
static field to allow the electrons to move towards the upper energy valleys, the Gunn 
field Eg = 3.25 kV/cm, we find a positive effect of the field fluctuations. In fact, despite 
the error bars are large, our findings show that the addition of correlated fluctuations, 
characterized by IO^^tq < td < tq, can increase the value of the spin depolarization 
length L up to the 20% of Lq. This effect is maximum for td/tq ^ 1. For the reasons 
discussed above, even in the high field case, both for very high and very low values of 
noise correlation time td, the value of L approaches Lq. The presence of a positive 
effect of noise can be ascribed to the reduction of the electron occupation percentage 
in L-valleys, shown in panel (b) of Figure [HI This finding, which can be considered 
as a further example of Noise Enhanced Stability (NES) [131 EHl |36], leads a greater 
number of electrons to experience a spin-orbit coupling in F-valley at least of one order 
of magnitude weaker than that present in L-valleys [26] , causing a decrease of efficacy 
of the DP dephasing mechanism. 
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Figure 4. Ratio between ttie spin depolarization length L in the presence of noise 
and Lq, obtained in absence of noise, as a function of the ratio between the noise 
amplitude D^/^ and Fq. The ee-points are obtained through full calculations including 
the electron-electron scattering mechanism; the no ee-points are calculated without it. 
(a) Td = O.lro and Fq = I kV/cm, (b) td = tq and Fq = 6 kV/cm. 



Although for many decades it has been believed that the electron-electron (e-e) 
Coulomb scattering does not contribute to the spin relaxation/dephasing process, in 
the presence of inhomogeneous broadening in spin precession, each kind of scattering, 
including spin-conserving scattering, can cause irreversible spin dephasing. Moreover, 
the Coulomb interaction is fundamental to obtain the correct distribution of the electrons 
in the k-space [3l |32]. Recently, it has been shown that e-e scattering inclusion leads 
to a strong increase of spin relaxation time in bulk semiconductors, in a wide range of 
electric field amplitude, lattice temperature and doping density [29] . 

In order to quantify the effect of the Coulomb interaction on the spin depolarization 
process in the presence of external noise, we show a comparison between the results 
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obtained with the e-e scattering mechanism and without it, analyzing the behavior of 
L/Lq^ as a function of the ratio between the noise standard deviation D^/^ and the 
deterministic value of the driving field Fq (see Figure H]). We show the comparison 
at the values of noise correlation time corresponding with the maximum of the noise 
effect observed by neglecting the Coulomb interaction. Each panel of Figure H] shows the 
"ee-points", obtained through full calculations including the electron-electron scattering 
mechanism, and the "no ee-points", calculated without the electron-electron interaction. 
In panel (a): td = O.Itq and Fq = 1 kV/cm, with Lq = 127 fim and tq = 0.49 ns; in 
panel (b): T£, = Tq and Fq = 6 kV/cm, with Lq = 666 nm and Tq = 2.42 ps; the 
new values of tq and Lq have been obtained by including the Coulomb scattering. For 
To = O.Ito, the reduction of the spin depolarization length, caused by fluctuations, 
is slightly affected from the inclusion of electron-electron scattering mechanism up to 
values of D^^"^ lower than 40% of the value of Fq. For values of D^^'^/Fq greater than 
0.4, the inclusion of the e-e scattering leads to a longer spin dephasing length, i.e. the 
electron-electron scattering mechanism reduces the negative effect of noise. This effect 
of the Coulomb interaction could be ascribed to the frequent momentum redistribution 
experienced from the electrons ensemble |3]. Unfortunately, the Coulomb interaction 
inclusion seems to randomize the system also for field amplitude greater than the Gunn 
field. Under high electric field conditions, in fact, up to values of D^^'^ lower than 33% 
of the value of Fq, the quantity L/Lq is almost not influenced by the inclusion of the 
e-e mechanism. By increasing the noise amplitude, a slight noise-induced positive effect 
on spin relaxation length is found. In this case, the addition of a source of correlated 
fluctuations, having correlation time comparable with the spin lifetime, enhances the 
value of the spin depolarization length L of only about 10-15%. 

In this work, for the first time, we have investigated the noise influence on the 
electron spin relaxation process in lightly ra-doped GaAs semiconductor bulks by also 
including the electron-electron interaction. The findings show that a fluctuating electric 
field, obtained by adding a correlated source of noise to a static field, can modify the 
spin depolarization length. For electric fields lower than the Gunn field and values of the 
noise correlation time T£, ~ Tq, the spin lifetime obtained in absence of noise, a reduction 
of the spin depolarization length up to 15% has been observed, strongly dependent on 
the noise intensity. This behavior can be explained by the different effective electric 
field experienced by the electron ensemble, within a time window comparable with tq. 
On the contrary, in the high electric field regime, for td=tq, we find an enhancement 
of the spin relaxation length up to 15%. This positive effect can be explained by the 
decrease of the occupation of the L-valleys, where the strength of spin-orbit coupling 
felt by electrons is at least one order of magnitude greater than that present in F-valley 
and represents an example of NES in spin depolarization process. 

To conclude, our preliminary results show that the presence of fluctuations in 
applied voltages changes the maintenance of long spin depolarization lengths in a way 
strongly dependent on both the strength of the applied electric field and the noise 
correlation time. Further studies are needed to learn more about the relationship 
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between the semiconductor characteristic time scales, the noise correlation time and 
the e-e interaction, in order to find the most favorable conditions for the manipulation 
of electron spins. 
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